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ABSTRACT 
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Context. Wide-range spectral coverage of blazar-type active galactic nuclei is of paramount importance for understanding the particle acceleration 
mechanisms assumed to take place in their jets. The Major Atmospheric Gamma Imaging Cerenkov (MAGIC) telescope participated in three 
multiwavelength (MWL) campaigns, observing the blazar Markarian (Mkn) 421 during the nights of April 28 and 29, 2006, and June 14, 2006. 
Aims. We analyzed the corresponding MAGIC very-high energy observations during 9 nights from April 22 to 30, 2006 and on June 14, 2006. We 
I , inferred light curves with sub-day resolution and night-by-night energy spectra. 
t-H ' Methods. MAGIC detects y-rays by observing extended air showers in the atmosphere. The obtained air-shower images were analyzed using the 
(^) standard MAGIC analysis chain. 

, Results. A strong y-ray signal was detected from Mkn 421 on all observation nights. The flux (E > 250 GeV) varied on night-by-night basis 
t-H ! between (0.92 ±0.11) lO^cm-V 1 (0.57 Crab units) and (3.21 ± 0. 15) 10- 10 cnr V 1 (2.0 Crab units) in April 2006. There is a clear indication for 
intra-night variability with a doubling time of 36 ± 10 sm , minutes on the night of April 29, 2006, establishing once more rapid flux variability for 
this object. For all individual nights y-ray spectra could be inferred, with power-law indices ranging from 1.66 to 2.47. We did not find statistically 
significant correlations between the spectral index and the flux state for individual nights. During the June 2006 campaign, a flux substantially 
lower than the one measured by the Whipple 10-m telescope four days later was found. Using a log-parabolic power law fit we deduced for some 
■ data sets the location of the spectral peak in the very-high energy regime. Our results confirm the indications of rising peak energy with increasing 
flux, as expected in leptonic acceleration models. 

Key words. Gamma rays: galaxies - BL Lacertae objects: individual (Mkn 421) - Radiation mechanisms: non-thermal 
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1. Introduction 

The active galactic nucleus (AGN) Markarian (Mkn) 421 was the 
first extragalactic source detected in the TeV energy range, using 
imaging atmospheric Cer enkov telescopes (IACTs: [Punch et all 
1992j iPetrv et all Il996l) . With a redshift of z = 0.030 it is 
the closest known and, along with Mkn 501, the best-studied 



TeV y-ray emitting blazarQ So far, flux varia t ions b y more than 
one order of magnitude (e.g., iFossati etail 120081). and occa- 
siona l flux doubling times as s hort as 15 min (iGai dos et al., 
1996; lAharonian et all l2002t ISchweizer. Wagner & Lorena, 
2008) have been observed. Variations in the hardness of 
the TeV y-ray spectrum during flares were reported by sev- 
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eral groups (e.g. Krennrich et all l2002t lAharonian et al 
Foss atiet all l2008). Simultaneous observations in the X-ray and 
very-high energy (VHE; E ^ 100 Ge V) bands show strong ev- 
idence for correla t ed flux variability (iKrawczvnski et ail [200 1 ; 
iBIazeiowski et all 120051; iFossati et all 120081) . With a long his- 
tory of observations, Mkn421 is an idea l candidate for long-term 
and statistical studies of its emission dTluczvkont et all l2007t 
iGoebel et alll2008atlHsu et al.Ll2009t) . 

Mkn 42 1 has been detected and studied at basically all wave- 
lengths of the electromagnetic spectrum from radio waves up to 
VHE y-rays. Its wide-range spectral energy distribution (SED) 
shows the typical double-peak structure of AGN. Mkn 421 is a 
so-called blazar. These constitute a rare subclass of AGNs with 
beamed emission closely aligned to our line of sight. In blazars, 
the low-energy peak at keV energies is thought to arise domi- 
nantly from synchrotron emission of electrons, while the origin 
of the high-energy (GeV-TeV) bump is still debated. The SED is 
commonly interpreted as being due to the beamed, non-thermal 
emission of synchrotron and inverse-Compton radiation from 
ultrarelativistic electrons. These are assumed to be accelerated 
by shocks moving along the jets at relativistic bulk speed. 
For most of the observations, the SED can be reasonably well 
described by homoge neous one-zone sy n chrotr o n-self-Compton 
(SSC ) models (e.g. iMarscher & GeaA [19851; iMaraschi et ail 
1 19921: ICo stamante & Ghise liiniL l2002l). Hadronic models 
(Mannhei m et all [1996; Miic ke et all 12003), however, can also 
explain the observed features. A way to distinguish between 
the different emission models is to determine the positions, 
evolution and possible correlations (see, e.g., Wagner, 2008b, 
for a review) of both peaks in the SED, using simultaneous, 
time-resolved data covering a broad energy range, e.g., as 
obtained in multiwavelength (MWL) observational campaigns. 

In this Paper we present results from Major Atmospheric 
Gamma-ray Imaging Cerenkov (MAGIC) telescope VHE y-ray 
observations of Mkn 421 during eight nights from April 22 to 
30, 2006, and on June 14, 2006. For most of the days, op- 
tical 7?-band observations were conducted with the KVA tele- 
scope. Simultaneous observations were performed by Suzaku 
(Mits uda et all l2007h and H.E.S.S., as well as by XMM-Newton 
Oansen et all 1200 lh on April 28 and 29, 2006, respectively. 
During both nights, we carried out particularly long, uninter- 
rupted observations in the VHE energy band of a; 3 hours du- 
ration each. An onset of activity in the X-ray band triggered 
an INTEGRAL-led target-of-opportunity (ToO) campaign, which 
took place from J une 14 - 25, 2006 for a total of 829 ks 
dLichti et all l2008h . Within this campaign, MAGIC observed 
Mkn 421 at rather high zenith angles from 43 to 52 degrees in 
parallel with INTEGRAL on June 14, 2006. 

In the following sections, we describe the data sets and the 
analysis applied to the VHE y-ray data, the determination of 
spectra for all observation nights, and put the results into per- 
spective with other VHE y-ray observations of Mkn 421. The 
interpretation of the se data in a MWL context is presented in 
lAcciari et all ((2009) and subsequent papers. 

VHE y-ray observations in April and J une 2006 have als o 
been carried out by the Whipple telescope dHoran et all 12 009), 
by th e VERITAS dFeganl [2008), and TACTIC dYadav et all 
2007) collaborations, although not simultaneously with our ob- 
servations. 



2. The MAGIC telescope 

The VHE y-ray observations were conducted with the MAGIC 
telescope located on the Canary island La Palma (2200 m above 



sea level, 28°45'N, 17°54'W). At the time of our observations 
in 2006, MAGIC was a single-dish 17-m instrument for 
the detection of atmospheric air showers induced by y-rays. Its 
hexagonally-shaped camera with a field of view (FOV) of « 3.5° 
mean diameter comprises 576 high-sensitivity photomultiplier 
tubes (PMTs): 180 pixels of 0.2° surround the inner section of 
the camera of 394 pixels of 0.1° (= 2.2° FOV). The trigger is 
formed by a coincidence of > 4 neighboring pixels. Presently the 
acces sible trigger energ y range (using the MAGIC standard trig- 
ger; iMeucci etall2 007) spans from 50 - 60 GeV (at small zenith 
angles) up to tens of TeV Further details, tele scope parame- 
ters, and performance information can be found inlBaixeras et all 
(l2004h : ICortina et alJ d2005l) : I Albert et all d2008al) . 



3. Observations and data analysis 

The observations were carried ou t during dark nigh ts, employ- 
ing the so-called wobble mode (iDaum et all \l99% . in which 
two opposite sky directions, each 0.4° off the source, are tracked 
alternatingly for 20 minutes each. The on-source data are de- 
fined by calculating image parameters with respect to the source 
position, whereas background control ("off") data are obtained 
from the same data set, but with image parameters calculated 
with respect to three positions, arranged symmetrically to the 
on-source region with respect to the camera center. The simulta- 
neous measurement of signal and background makes additional 
background control data unnecessary. In order to avoid an un- 
wanted contribution from source y-events in the off sample, and 
to guarantee the statistical independence between the on and the 
off samples in the signal region, events included in the signal re- 
gion of the on sample were excluded from the off sample and 
vice versa. 

The data we re analyzed following the standard MAGIC anal - 
ysis procedure dBretz & Wagneii 120031; iBretz & Dorneit 12008). 
After calibration ( Alber t et al. , 2008c) and extracting the signal 
at the pulse maximum using a spline method, the air-shower im- 
ages were cleaned of noise from night-sky background light by 
applying a three-stage image cleaning. The first stage requires 
a minimum number of 6 photoelectrons in the core pixels and 
3 photoelectr ons in the boundary pixels of the images (see, e.g. 
FegaBtHll. These tail cuts are scaled according to the larger 
size of the outer pixels of the MAGIC camera. Only pixels with 
at least two adjacent pixels with a signal arrival time difference 
lower than 1 .75 ns survive the second cleaning stage. The third 
stage repeats the cleaning of the second stage, but requires only 
one adjacent pixel within the 1.75 ns time window. 

The data were filtered by rejecting trivial background events, 
such as accidental noise triggers, triggers from nearby muons, or 
data taken during adverse atmospheric conditions (e.g., low at- 
mospheric transmission). 12.7 hours out of the total 15.0 hours' 
worth of data survived the latter quality selection and were used 
for further analysis. 

We calculated image parameters dHillasL 1 19851) such as 
WIDTH, LENGTH, SIZE, CONC, M3LONG (the third mo- 
ment of the light distribution along the major image axis), 
and LEAKAGE (the fraction of light contained in the outer- 
most ring of camera pixels) for the surviving events. For the 
y/hadron separation, a SIZE-depend ent parabolic cut in AREA = 
WIDTH x LENGTH x n was used dRiegel et alll2005l) . The cut 
parameters for the assessment of the detection significance were 
optimized on Mkn 421 data from close-by days. For the data of 
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June 14, 2006 at rather large zenith angles, data of Mkn 501 from 
October 2006 were used to determine the optimal cuts. Any sig- 
nifica nce in this work was calculated using Eq. 17 of iLi & Ma 
(119831) with a = 1/3. 

The primary y-ray energies were reconstructed from the 
image param e ters us ing a Random Forest regression method 
(lAlbert et all l2008bl and references therein) traine d with 
Monte-Carlo simulate d events (MCs; iKnapp & HeckL 12004 
Maiumda r et all 120051) . The MC sample is characterized by a 
power-law spectrum between 10 GeV and 30 TeV with a dif- 
ferential spectral photon index of a — -2.6, and a point-spread 
function resembling the experimental one. The events were se- 
lected to cover the same zenith distance range as the data. For the 
spectrum calculation, the area cut parameters were optimized to 
yield a constant MC cut efficiency of 90% over the whole energy 
range, increasing the y-ray event statistics at the threshold. 

The Mkn 421 observations presented here are among the first 
data t aken by MAGIC afte r major hardware updates in April 
2006 (iGoebel et all l2008bl) . which required us to thoroughly 
examine the data. Despite the hardware changes, the MAGIC 
subsystems performed as expected with the exception of an un- 
stable trigger behavior for some PMTs, leading to a significant 
loss of events in one of the six sectors of the camera. In order 
to proceed with the data analysis with serenity and to estimate 
the effect caused by this inhomogeneity, a simple procedure was 
applied to the data: The expected number of events, as a func- 
tion of energy, for the affected sector was estimated as the mean 
of the number of events in the other five sectors of the camera. 
(A homogeneous distribution of events through the six sectors 
is expected for normal conditions). The difference between the 
expected and actually measured events was computed using the 
whole data sample in order to have sufficient statistics. We found 
a decrease of the differential photon flux of 5.7% between 250 
and 400 GeV, 4.6% between 400 and 650 GeV, 2.2% between 
650 and 1050 GeV and < 1% for higher energies for the April 
2006 data. Due to the higher zenith distance and energy thresh- 
old, the method was adapted for June 14, 2006 and yielded a de- 
crease of 5.2% between 450 and 670 GeV and 2.6% for higher 
energies. However the above mentioned effect is just an average 
one, with estimated flux errors of up to 6.6% showing up for 
individual nights. 

To mitigate the effect of the inhomogeneity, instead of an (al- 
ready increased) energy threshold of 250 GeV, higher thresholds 
of 350 or 450 GeV were applied for some observation nights. 
In this way we made sure that the estimated systematic error re- 
mains within reasonable limits. 

For the calculation of the individual light curves as well as 
for the overall April 2006 lightcurve, the flux between 250 GeV 
and 350 GeV was extrapolated for the nights with higher thresh- 
old. We assumed a power-law behavior in this energy range, with 
the spectral index determined for the first three energy bins of 
the whole April dataset (i.e., a = -2.08). The flux normalization 
for each night has been determined at 500 GeV by a fit to the 
first three differential spectral points, an energy range which is 
reliable for all affected nights. 

Tab.Q]summarizes the analyzed data sets. The statistical sig- 
nificance of any detection is assessed by applying a cut in 6 2 , 
where 9 is the angular distance between the expected source po- 
sition and the reconstructed y-ray arrival direction. The arrival 
directions of the showers in e quatorial coordinates were calcu- 
lated using the DISP method (Fom in et al l ll994ULessardetall 
1200 ll) . We replaced the constant coefficient £ in the parameteri- 



Table 1. Some characteristic parameters of the different data 
sets of the campaign. 



Night 


Observation Window [MJD] 


h«. [h] 


ZAP] 


April 22, 2006 


53847.97679 


- 53848.01460 


0.76 


18 


-28 


April 24, 2006 


53849.96428 


- 53850.00669 


0.99 


16 


-28 


April 25, 2006 


53850.92813 


- 53850.99607 


1.53 


10 


-26 


April 26, 2006 


53851.92862 


- 53852.00383 


1.64 


10 


-29 


April 27, 2006 


53852.93474 


- 53853.00047 


1.42 


12 


-28 


April 28, 2006 


53853.88173 


- 53854.01394 


2.23 


10 


-32 


April 29, 2006 


53854.89514 


-53855.04119 


2.78 


9 


-41 


April 30, 2006 


53855.97283 


- 53855.97906 


0.16 


23 


-24 


June 14, 2006 


53900.91979 


- 53900.95532 


0.80 


43 


-52 



Notes. f cff denotes the effective observation time. ZA gives the zenith 
angle range of the observations. 



zation of DISP in the original approach by a term which is de- 
pendent on LEAKAGE, SIZE, and SLOPE, 

^ = ^o+^i SLOPE + ft LEAKAGE + k& (log 10 SIZE - &) 2 , (1) 

k — for log I0 SIZE < £4 and k — 1 for log I0 SIZE > £4. The co- 
efficients were determined using simulated data. The parameter 
SLOPE is a measure for the longitudinal arrival time evolution 
of the shower i n the camera plan e similar to the time parameter 
GRADIENT in lAliu et al.ld2009l) . Instead of defining the param- 
eter from a fit to the arrival time distribution, however, SLOPE is 
determined as an analytical solution of the fit. Note that this new 
parametrization makes DISP and therefore 9 2 source dependent. 

All stated errors are statistical errors only; we estimate our 
systematic errors to be 16% for the energy scale, 11% for ab- 
solute fluxes and flux norm alizations, and 0.2 for the spectral 
slopes (lAlbert et all 1200 8a), not including the additional sys- 
tematic flux errors mentioned above. 

A second, independent analysis of the data yielded compati- 
ble results to those presented here. 

4. Results 

4. 1 . Results for April 22 - 30, 2006 

MAGIC observed Mkn421 from MJD 53847 to MJD 53855. 
During the observations, two MWL campaigns were carried out 
simultaneously with Suzaku and with XMM-Newton on MJD 
53854 and MJD 53855, respectively. Mkn 421 was also observed 
as par t of the monitorin g program of the Whipple 10-m telescope 
(see lHoran etafll2009h . albeit about 3.5 hours after the MAGIC 
observations stopped, due to the different longitudes of the two 
instruments. 

A strong y-ray signal from the source was detected in all 
eight observation nights. In total, 3165 excess events were 
recorded over a background of 693 events for energies > 
250 GeV, yielding an overall significance of 64.8<x. Mkn 421 
exhibited an average flux of f>250GeV = (1.48 + 0.03) • 
10" 10 cm~ 2 s~ 1 . When compared to earlier observations (see, 
e.g. lAlbert et all l2007at [Tluczvkont et all l2007t fcoebel et all 
2008at lSteele et alll 2008). our observations indicate an elevated 
flux state of Mkn 421. We found high flux states in the nights of 
MJD 53850, F>250GeV = (2-32 + 0.13) ■ lO-^crrrV 1 , MJD 
53853, F >25 oGeV = (3-21 + 0.15) ■ lCr^cnrV 1 , and MJD 
53856, F >25 oGeV = (2.39 + 0.33) • 10" 10 cnrV 1 (Fig. H). In 
the remaining nights (we assumed nights with fluxes below 
1.6 • 10~ 10 cm~ 2 s _1 as non-flare nights), Mkn 421 exhibited a 
low-flux average of F >25 oGev = (1-09 + 0.03) ■ lCT^cm-V 1 . 
The analysis results on a night-by-night basis are summarized 
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Fig.l. VHE ( E > 25 GeV) light curve for Mkn421 observations in April 2006. The dotted line represents the Crab nebula flux 
(Albe rt et all l2008ah . whereas the individual dashed lines show the result of a fit to the time bins (average nightly flux) of the 
corresponding nights. 
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Fig. 2. VHE (E > 250 GeV) light curve for Mkn421 observa- 
tions in April 2006. The data points represent average nightly 
fluxes. The observation windows of the Suzaku (MJD 53853.28- 
53854.27) and XMM-Newton (MJD 53854.87-53855.35) MWL 
campaigns are marked by the gray-shaded areas. A "mean 
low flux" (solid line) was averaged over all data points below 
1.6 • 10 _10 cirr 2 s _1 , i.e., excluding those data points marked by 
thin open circles. T he dashed line gives the Crab nebula flux 
(lAlbertet all 12008 a) for comparison. 



in Tab. [2] and include the nightly numbers for excess and back- 
ground events, significances, and average integral fluxes above 
250 GeV (where the nights with an energy cut of 350 GeV where 
extrapolated down to 250 GeV, see Sect. [3]for details). The re- 
sults of a spectral fit based on a simple power law (PL) of the 
form 



dF 
dE 



fo ■ lO^TeV-WV 1 



(2) 



are also shown. 

The energy thresholds of the individual observations are also 
given in Tab. [2] As the analysis threshold is always lower than 
the applied energy cut, the latter one defines the energy threshold 
value. 

The strong y-ray signal allowed to infer light curves with 
a resolution below one hour for all of the observation nights, 
which are shown in Fig. [T] (see Tab. [3] for the light curve data). 
Most light curves are compatible with a constant flux during 



the nightly observation time (see Tab. [2] for all constant-fit 
values), while on MJD 53855 a clear intra-night variability is 
apparent. A fit with a constant function yields an unacceptable 
*red = 41/7 (P « 8 • 10~ 5 %) for this night, and the data suggest 
a flux halving time of 36 + 10 sfa/ minutes. Note that this interest- 
ing observation window has als o been covered by XM M-Newton 
observations in the X-ray band (lAcciari et all 2009). 

4.2. Results for June 14, 2006 

An onset of activity to « 2 times the average quiescent-flux 
level of Mkn421 was measured in April 2006 by the RXTE 
all-sky monitor (ASM) instrument. It triggered an INTEGRAL 
ToO campai g n from June 14, 2006 to 25 for a total of 829 ks 
dLichti et all l2008h . This > 30 mCrab flux remained until 
September 2006. During the 9-day campaign, Mkn421 was tar- 
geted by various instruments in the radio, optical, X- r ay an d 
VHE wavebands. Results are reported in iLichti et aD {2008). 
On June 14, 2006, MAGIC observed Mkn421 at rather high 
zenith angles in parallel with the OMC, JEM-X, and IBIS mea- 
surements aboard INTEGRAL. Further VHE coverage was pro- 
vided by the Whip ple 10-m telescope on June 18/19/21, 2006 
dLichti et all [2008b . 

The MAGIC observations on June 14, 2006 lasted for « 50 
minutes. The high zenith angles of 43 to 52 degrees of this obser- 
vations and the previously mentioned inhomogeneities result in 
an energy threshold of £Vesh. - 450 GeV. In spite of the overall 
rather difficult observatio nal circumstances caused by the hi gh 
zenith angle observations dTonellol 120061; lAlbert et"all l2006b . a 
firm detection on the 7.5-cr significance level was achieved. 

The corresponding differential energy spectrum is shown in 
Fig. [3] Between 450 GeV and 2.2 TeV, it can be described by a 
simple power-law of the form 



— = (1.68 + 0.32)- 10 -11 TeV _1 cnTV 
dE 



ll.OTeV/ 



-2.38±0.44 



(3) 



For comparison we also show the spectral points reported 
by the Whipple 10-m telescope averaged over the nights of 
June 18/19/21, 2006. Generally, there might be systematic dif- 
ferences between the Whipple and MAGIC measurements. It 
could, however, be shown that such inter-instrument systematic 
effects are rather sma ll and under contro l, e.g. those between 
MAGIC and H.E.S.S. dMazin et alll2005h . Particularly the Crab 
nebula spectra measured by Whipple and MAGIC agree quite 
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Table 2. Analysis results. 



Observation Night 


excess 


Nbackgr. 


S 


E cu , [GeV] 


F(E > E mm ) 


i 

A red.const 




fo 




a 




2 

Ared.PL 


April 22, 2006 


100 


29 


10.9cr 


350 


0.92 = 


t0.ll 


1.3/2 


0.98 


+ 0.15 


2.05 


± 


0.21 


2.1/2 


April 24, 2006 


419 


69 


25.0a- 


250 


2.32 = 


= 0.13 


2.7/2 


2.45 


+ 0.14 


2.25 


± 


0.09 


2.0/3 


April 25, 2006 


342 


83 


20.8a- 


250 


1.34 = 


t0.09 


1.7/2 


1.43 


±0.09 


2.26 


± 


0.12 


0.24/3 


April 26, 2006 


225 


62 


16.4a- 


350 


1.08 = 


t0.09 


1.3/4 


1.21 


+ 0.11 


2.35 


± 


0.17 


0.41/2 


April 27, 2006 


615 


56 


33.5a- 


350 


3.21 = 


1=0.15 


1.9/4 


3.37 


±0.18 


2.07 


± 


0.07 


4.8/4 


April 28, 2006 


311 


75 


19.9a- 


350 


1.14 = 


t0.08 


4.3/8 


1.32 


±0.10 


2.47 


± 


0.14 


0.65/2 


April 29, 2006 


514 


169 


23.7a- 


250 


1.04 = 


= 0.06 


41/7 


1.14 


±0.06 


2.28 


± 


0.09 


2.0/4 


April 30, 2006 


69 


11 


10.3a- 


250 


2.39 = 


t0.33 




2.16 


±0.34 


1.66 


± 


0.20 


1.4/1 


June 14, 2006 


95 


87 


7.5a- 


450 


0.34 = 


t0.06 


2.4/1 


0.168 


± 0.032 


2.38 


± 


0.44 


1.5/2 



Notes. Number of excess (iV cxccss ) and background (A'backgr.) events, resulting significances 5, lower cuts in event energy, integral fluxes F above 
E mm = 250 GeV for the April 2006 data and E mm = 450 GeV for the June 14, 2006 data (in units of 10~ 10 cnrV), fit quality of a constant-flux 
fit to the individual observation nights (see Fig. [T}, and power-law fit results for the differential energy spectra of dF/dE = f ■ (E/E y with 
E = 0.5 TeV for the April 2006 data and E = 1.0 TeV for the June 14, 2006 data, respectively; f in units of lO^TeV^cnrV 1 . 



Table 3. Light curve data. 



Observation 

[MJD] [10- 1C 
2006/04/22 



53847.98307 UW: 

53847.99775 0.70 : 
53848.00867 0.99 : 

2006/04/24 

97136 



0.23 
0.17 



~33W 

53849 
53850 



236": 

98618 2.04: 
00033 2.37 : 
— 2006/04/25 



U72T" 
0.21 
0.22 



53850. 
53850. 
53850. 



93996 1.24: 
96431 1.49: 
98652 1.26: 
— 2006/04/26 



0.15 
0.15 
0.15 



53851.93677 


0.97 


±0.20 


53851.95255 


1.04 


±0.21 


53851.96726 


1.25 


±0.20 


53851.98190 


1.00 


±0.18 


53851.99680 


1.13 


±0.18 



2006/04/27 



53852.94098 


3.01 = 


t 0.33 


53852.95502 


3.19 = 


t0.38 


53852.96823 


3.05 = 


t0.31 


53852.98159 


3.57 = 


t0.32 


53852.99406 


3.17 = 


t0.28 



53853.88754 


1.36 


±0.25 


53853.89880 


0.80 


±0.27 


53853.92893 


1.22 


±0.24 


53853.93984 


1.09 


±0.25 


53853.95457 


1.18 


±0.24 


53853.96887 


1.22 


±0.27 


53853.98040 


0.95 


±0.19 


53853.99316 


1.32 


±0.21 


53854.00687 


1.12 


±0.18 



2006/04/29 



53854.90199 


2.07 = 


tO.23 


53854.91620 


1.42 = 


t0.19 


53854.95206 


0.86 = 


t0.16 


53854.96625 


1.11 = 


= 0.17 


53854.97974 


1.27 = 


t 0.18 


53854.99354 


0.83 = 


t 0.15 


53855.00847 


0.80 = 


= 0.14 


53855.02879 


0.69 = 


t0.ll 



53855.97595 



"33900" 

53900. 



2.39 ± 0.33 
2006/06/14 



9T79T" 
94585 



H4T: 
0.26: 



0.08 



10" 



,.10" 




•This Work, June 14, 2006 
o Lichti et al. (2008), 
June 18/19/21, 2006 



600 



1000 
Energy [GeV] 



2000 



Fig. 3. Differential photon spectrum for Mkn 421 for the obser- 
vation night of June 14, 2006 (black data points). A power-law 
fit to the spectrum results in a spectral slope of a = -2.38 + 0.44 
(See Tab. [2] for the fit results). Also shown are spectral points 
measured with the Whipple 10-m telescope (Lichti et al., 2008) 
during June 18-21,2006. 



well (lAlbert et al.L l2008al) . The Mkn 421 flux measured by the 
Whipple 10-m telescope four days after the MAGIC observation 
is substantially higher than our measurements (Fig. [3]), point- 
ing to a clear evolution of the source emission level within the 
INTEGRAL campaign. 



5. Discussion 

In leptonic acceleration models, e.g., SSC models, a shift of the 
high-energy peak (attributed to Inverse Compton radiation) in 
the spectral energy distribution towards higher energies with an 
increasing flux level is expected. In the VHE domain, such a shift 
can be traced by spectral hardening. Variations in the hardness 
of the TeV y-ray spectrum during flares were reported by sev- 
eral groups (e.g.. [ Krennrich et al.Ll2~002t lAharonian et all [2005 ; 
iFossati et all 120081) . We tested for a correlation of the spectral 
hardness with the flux level of the de-absorbed spectrum (i.e. af- 
ter removing any attenuation effects caused by the Extragalactic 
Background Light fEBLl. cf.lNikishovll 19621: [Gould & Schrederl 
Il966t lHauser&Dwekl 120011) in our data (Fig. gji, but found 
that the correlation neither can be described by a constant fit 
(X^ ed = 17/8, P ^ 3%) nor by a linear dependence of spectral 
hardness and flux level (xLi = 11/7, P « 12%), giving no clear 
preference for either. Although clear flux variations are present 
in the data set, the overall dynamical range of 3.9 in flux might 
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Fig. 4. Spectral index vs. flux at 0.5 TeV deduced from a sim- 
ple power-law fit after EBL de-absorption for Mkn 42 1 . The x^ ed 
for a constant fit (spectral index uncorrelated with flux level; 
solid line) amounts to 17/8 (P ~ 3%), while a linear correla- 
tion (dashed line) has a^ 2 , = 11/7, equal to P m 12%. 



be too small to see a significant spectral hardening with increas- 
ing flux. 

The individual night-by-night spectra during the campaign 
in April 2006 are shown in Fig. [5] All spectral data points are 
summarized in Tab. [4] For the nights of April 22, 26, and 29, 
2006, there seems to be evidence for a resolve d peak, but a like- 
lihood ratio test (e.g. jMazin & Gqebell l2007h yields significant 
curvature only for April 27, 20060 We used a logarithmic cur- 
vature term, corresponding to a paraboli c power-law (log-P) i n 
a \og(E 2 dF/dE) vs. logE representation dMassaro et al.L 12004). 
and a power-law with exponential cutoff (PL+C) of the form 



dF -f 
dE~ f ° 



and 



10 TeV 



1 -2 - 

cm s 



E \ -^+,eiog 10 ^j 
E~ ) 



(4) 



g =/0 . 1 0-.TeV-e m -v(|)" O e X p(^), 

respectively. The likelihood ratio test results in a clear preference 
towards a log-P or a PL+C compared to a simple power-law with 
a probability of s* 96% for both of them. The ^ 2 ed values for 
PL, log-P, and PL+C fits on the individual night-by-night spectra 
in Fig. [5] are given in Tab. [5] Also the high statistics data sets 
defined by combining all data from April, all data from the five 
low-state nights and all data from the three high-state nights, 
clearly showed evidence for a parabolic or cutoff shape of the 
spectra. The results of the fits and the probability of a likelihood 
ratio test are given in Tab. [6] For all these nights our data did 
not allow to prefer one model over the other. The fact that all of 
the high statistics data sets show a curved spectral shape is an 
indication of this feature being always visible for Mkn 421 and 
hence source intrinsic. 

The curved power laws enable to locate a peak in the de- 
absorbed spectrum at £p ea k = EqIQ^"^^ for the log-P and at 
£peak = (2 - a)E cut if a < 2 for the PL+C fit. For simplicity 



Table 4. Energy spectra for all observation nights under study 
after EBL de-absorption. 



E bounds 
[GeV] 



Flux 
[TeV cm -2 s 



2006/04/22 


350 


554 (2.39 + 0.56) 


10~ 




554 


877 (3.67 ± 0.75) 


io- 




877 


1389 (3.29 + 0.79) 


io- 




1389 


2200 (2.51 +0.93) 


io- 




2006/04/24 


250 


435 (7.00 ± 0.64) 


io- 




435 


758 (7.69 ± 0.76) 


io- 




758 


1320 (6.12 + 0.86) 


io- 




1320 


2297 (6.54 ±1.33) 


io- 




2297 


4000 (4.80 ±1.74) 


io- 




2006/04/25 


250 


416 (4.36 ±0.43) 


io- 




416 


693 (3.95 ± 0.48) 


io- 




693 


1154 (3.86 ±0.62) 


io- 




1154 


1922 (3.47 ± 0.87) 


io- 




1922 


3200 (3.85 ± 1.31) 


10" 




2006/04/26 


350 


572 (3.41 ±0.41) 


io- 




572 


935 (3.46 ± 0.49) 


io- 




935 


1529 (2.90 ± 0.65) 


io- 




1529 


2500 (2.38 ± 0.75) 


io- 




2006/04/27 


350 


549 (8.83 ±0.66) 


io- 


1 


549 


860 (1.07 ±0.09) 


io- 





860 


1349 (1.09 ±0.12) 


io- 





1349 


2115 (1.17 + 0.17) 


io- 





2115 


3317 (1.03 ±0.26) 


io- 





3317 


5200 (6.89 ± 2.76) 


io- 


1 


2006/04/28 


350 


635 (3.92 ± 0.34) 


io- 




635 


1153 (3.01 ±0.39) 


io- 




1153 


2093 (2.96 ± 0.57) 


io- 




2093 


3800 (2.06 ± 0.79) 


10" 




2006/04/29 


250 


387 (3.37 ±0.32) 


io- 




387 


600 (3.24 ± 0.32) 


io- 




600 


929 (3.29 ± 0.39) 


io- 




929 


1438 (3.33 ±0.45) 


io- 




1438 


2228 (2.31 ±0.53) 


io- 




2228 


3450 (2.08 ±0.81) 


io- 




2006/04/30 


250 


572 (4.87 ± 1.09) 


io- 


1 


572 


1310 (1.01 ±0.21) 


io- 





1310 


3000 (1.01 ±0.33) 


io- 





2006/06/14 


450 


669 (2.76 ± 0.74) 


io- 


1 


669 


995 (1.61 ±0.67) 


io- 


1 


995 


1480 (2.54 ± 0.79) 


io- 


1 


1480 


2200 (1.80 ±0.85) 


io- 


1 



Notes. The two energy bounds specify the range in which the corre- 
sponding flux was measured. 

Table 5. ^- 2 ed values for the PL, log-P, and PL+C fits performed 
in Fig. 





22 


24 


25 26 


27 


28 29 


30 


PL 


2.2/2 


1.9/3 


0.21/3 0.47/2 


5.3/4 


0.59/2 2.3/4 


1.5/1 


log-P 


0.14/1 




0.041/1 


0.48/3 


1.1/3 




PL+C 


0.27/1 




0.076/1 


0.34/3 


0.87/3 





Notes. The columns represent days in April 2006. 



3 the respective log-P probabilities are 83%, 48%, 73%, and 96%. 
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Fig. 5. Differential energy spectra for Mkn421 for April 2006 before (gray points) and after (black points) correcting for EBL 
absorption. For the apparently hard spectra on April 22, 26, 27, and 29, 2006, log-P (Eq. [4]i and PL+C (Eq. [5]) fits were performed 
(red solid and blue dashed curves, respectively). 



we determined £p ea k of the log-P by using the apex form of the 
parabola in a logarithmic representation: 

1Ogl0 cif = logl0 ^° +1 ° gl ° Q ''( l0gl(, (^/^f )) (6) 

which naturally yields both £p e ak and the flux at the peak, /o, re- 
spectively. Additionally, the spectral cutoff is naturally obtained 
from the PL+C fit as the fit parameter £ cllt . The results are shown 
in Tab. [6] The values of £p ea k as determined using the log-P and 
the PL+C were compatible with each other for the data sets av- 
eraging several nights and showed indications for an increase of 
the peak energy with rising flux level, as predicted if the VHE 



radiation were due to SSC mecha nisms. We compare our results 
with historical values taken from I Albert et alJ ll2007al) inFig.l6l 
Our data confirm the previously suggested correlation. 

The observation of a relation between flux (and thus, fiuence) 
and the position of the VHE peak in the SEP could be signa lling 
a relation simil ar to the one suggest ed bv lAmati et alJ d2002l) and 
observed by Sakamoto et al. (2006) for gamma-ray bursts. Since 
the TeV y-ray production is assumed to take place in a relativis- 
tic jet, and many of the same radiative processes are involved 
(on a larger scale, of course) it might be a similar (or related) 
mechanism at work on a different scale. A trend towards a re- 
lation between flux an d spectral index in the TeV energy range 
has also been noted by Wagner (2008a), studying 17 known TeV 
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Table 6. Special fit results. 

Data Set Used Fit fo c/> j3_ E cul [TeV] x^am Likelihood £ peak [TeV] 

PL 9.54 ±0.52 1.92 ±0.07 5.3/4 

log-P 9.35 ±0.55 1.54 ±0.19 0.59 ± 0.29 0.48/3 96% 1.2 ±0.7 

April 2/, IVUb i g.p a pex 11.5 ±0.9 0.26 ±0.17 0.48/3 96% 1.2 ±0.2 

PL+C 11.3 ±1.2 1.44 ±0.24 2.6 ± 1.3 0.34/3 96% 1.4 ± 1.0 

PL 4.53 ±0.07 2.07 ±0.04 16/5 

... . ., n log-P 4.75±0.12 1.89 ±0.06 0.39 ±0.11 1.2/4 99% 0.69±0.14 

All April uata i og .p ape x 4.84 ±0.16 0.41 ±0.11 1.2/4 99% 0.69 ± 0.06 

PL+C 5.36 ±0.31 1.77 ±0.09 3.6 ± 1.1 1.8/4 99% 0.80 ± 0.42 

PL 8.19 ±0.28 1.93 ±0.05 6.0/4 

W h <?t \r ht log " P 8.46 ±0.32 1.79 ± 0.09 0.29 ±0.15 2.0/3 94% 1.1 ±0.6 

Hign-Mate JNignts log . p apex 9,21 ±0.47 0.52 ±0.17 2.0/3 94% 1.1 ±0.3 

PL+C 9.02 ±0.64 1.75 ±0.12 6.1 ±4.0 3.1/3 91% 1.5 ± 1.2 

PL 3.39 ±0.10 2.17 ±0.05 6.6/4 

log-P 3.55 ±0.13 2.02 ±0.08 0.38 ±0.17 1.1/3 97% 0.48 ±0.12 

log-P apex 3.55 ±0.13 0.41 ±0.16 1.1/3 97% 0.48 ±0.12 

PL+C 4.15 ±0.40 1.85 ±0.15 2.9 ± 1.3 0.75/3 97% 0.45 ± 0.47 

Notes. Results of a power-law fit (Eq.[2}, a log-parabolic fit (Eq.f4]-(£/£o) 2 )> a log-parabolic fit in apex form (Eq.[6j and a power-law fit with an 
exponential cutoff ((Eq.[5j -(E/E Q ) 2 ) in E 2 dF/dE after EBL de-absorption for special data sets. f is given in units of 10~"TeV~ 1 cnT 2 s - '; a, a', ft, 
E cut and £peak are the fit parameters as stated in the text, and Likelihood denotes the probability of a likelihood ratio test. The on/off normalization 
factor is 1/3, E = 0.5 TeV. 
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Fig. 6. Derived peak position using the log-P (Eq. [6]i versus flux 
at 1 TeV f or the data se t s prese nted in Tab. [6] Historical data, 
taken from lAlbert et all d2007al) . are shown in gray. Our data 
confirm the indication of a correlation between the two parame- 
ters. 



blazars, and bv lTramacerd d200 9) in the X-ray band, after a deep 
spectral analysis of all Swift observations of Mrk 421 between 
April and July 2006. 

Although the peak energy measured on April 27, 2006 ex- 
ceeds that of the All April Data and Low-State data set, it is, 
despite having a higher flux, comparable with that derived for 
the High-State data set. This discrepancy in terms of the ex- 
pected behaviour in SSC models can be explained with the dif- 
ferent nature of the data sets: The April 27, 2006 data repre- 
sent a rather particular, 1 .4 h long episode of an individual flare 
event, whereas the High-State data set is an average of three in- 
dividual flares. Due to the sparse sampling, most probably each 
of these observations caught different epochs of the individual 
flare evolutions, during which the spectral shape can change 
considerably in terms of spectral index and curvature (see, e.g., 
iKatarzvhski et aUl2006l) . Hence the two data sets are not neces- 
sarily directly comparable. 

The values of the derived cutoff energies are also suggesting 
this behavior, showing, with the exception of April 27, 2006, 



10 



10 



£ 
o 

> 



10" 




SEDs from Literature 

• MAGIC (2004-2005) 
VERITAS (2008) 
HEGRA (1999-2000) 
HEGRA (2000-2001) 
HEGRA CT1 (2001) 
CAT (1998) 

• CAT (2000) 
Whipple (2000-2001), high 
Whipple (2000-2001), low 
HESS (2004) 



This Work, April 2006 

• High State 

• Low State 



_l_ 



_l_ 



_l_ 



i i 



J I 



10' 



10 J 

Energy [GeV] 



10 4 



Fig. 7. EBL de-abs orbed historical spectra of Mkn421 (see 
Albe rt et a 1., 2007 a, for references) along with selected spec- 
tra from the April 2006 campaign and the flare spectrum of 
Donnarumma et al .] (120091) . The solid line is the result of a fit 
using Eq. [4] No te that the historical d ata were deabsorbed us- 
ing the model of | Primack et al.l (12005), our data and those from 
iDonnaramma et al.1 (120091) with the model of iKneiske & Dolel 
(I2008h . 



an increase with rising flux, thus indicating a source-intrinsic 
rather than a cosmol ogical reason for the cuto ff feature. This is in 
accordance with the lKneiske & Dolel J2008) lower-limit model, 
predicting an EBL cutoff for Mkn421 at around 13 TeV. 

In Fig. |71 we compare "historical" spectra measured be- 
tween 1998 and 2005 with the low-state and high-state spec- 
tra derived from the observations reported here. It is obvious 
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that our low-state spectrum represents one of the lowest flux 
states ever measured in VHE for Mkn 42 1 , whereas the high state 
spectrum shows no exceptionally high flux level of this source. 
Both spectra are harder than historical spectra with compara- 
ble flux levels, in particular harder than the VERITAS spectrum 
dDonnarumma et al ., 2009), enabling one of the best measure- 
ments of the turnover of the SED in a low flux state. While a 
previou s observation yie l ded a rather flat spectrum in the VHE 
regime (Ahar onian et al.L |2002|) . we conclude that we measured 
a rather clear peak (flat structure in the SED). The low-state 
spectrum has a shape similar to the one measured by HEGRA 
CT1, although at an approximately three times lower flux level. 
The high-state spectral shape resembles the high-state Whipple 
spectrum, which in turn has an about three times higher flux. 
This tendency can also be seen in Fig. |6l which shows that the 
fluxes we derive are systematically lower than historical mea- 
surements for comparable peak energies. Within the SSC frame- 
work this difference in flux for comparable spectral shapes can 
be caused by, e.g., a lower number of electrons with the same 
energy distribution as in the high-flux case. 

In summary, we followed the evolution of a sequence of 
mild flares of the blazar Mkn 421 during one week from April 
22 to 30, 2006, peaking at F(E > 250 GeV) = (3.21 + 
0.15) 10 _10 cm _2 s _I (w 2.0 Crab units). The nocturnal obser- 
vations lasted at least for about one hour and allowed for the 
reconstruction of night-by-night spectra. During three observa- 
tion nights high fluxes were recorded, in which, however, no 
variability could be measured. In two of these nights, rather 
hard spectral indices were found, but this was also the case for 
the night with the lowest flux. During the night of April 29, 
2006, with a not particularly high flux of F(E > 250 GeV) = 
(1.04 + 0.06) lO ^cirrV 1 (« 0.65 Crab units), clear intra-night 
variability with a flux-doubling time of 36 + 10 s , fl , minutes was 
observed. 

According to a likelihood ratio test, the spectra of some data 
sets were better described by curved power laws than simple 
power laws, enabling us to calculate peak and cutoff energies 
in the VHE regime. The derived peak values are consistent with 
an evolution of the peak energy with the flux, as suggested by 
historical data. Indications of an intrinsic cutoff in the spectra of 
Mkn 421, as found in former observations, are confirmed by our 
results. 

During the INTEGRAL-triggeied MWL campaign in June 
2006 we observed Mkn 421 in one night at high zenith an- 
gles. Our measurements complement the three-night observa- 
tions conducted by the Whipple 10-m telescope four days later. 
Taking the MAGIC and Whipple results together, a variabil- 
ity of Mkn 421 also during the INTEGRAL observations is ev- 
ident. The energy coverage of the Whipple telescope spectrum 
(A £ « 600 GeV) was not sufficient to assess any spectral evo- 
lution by comparing it to the MAGIC spectrum (A E « 2 TeV). 

The determined fluxes and spectra will be further used for 
studies of the SED taking into account data taken at other photon 
energies in detailed MWL analyses (publications in preparation). 
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